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Abstract 

In this paper we study a simple example of a two-parameter space of renormalisa- 
tion group flows of defects in Virasoro minimal models. We use a combination of 
exact results, perturbation theory and the truncated conformal space approach 
to search for fixed points and investigate their nature. For the Ising model, we 
confirm the recent results of Fendley et al. In the case of central charge close 
to one, we find six fixed points, five of which we can identify in terms of known 
defects and one of which we conjecture is a new non-trivial conformal defect. We 
also include several new results on exact properties of perturbed defects and on 
the renormalisation group in the truncated conformal space approach. 
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1 Introduction 

By a defect in a two-dimensional conformal field theory we mean a line of inhomogeneity on 
the surface, where the expectation values of fields are allowed to be discontinuous or even 
singular. An example would be the continuum limit of a lattice model where the couplings 
are altered from their normal values along a line. 

A typical defect is not invariant under a scale transformation and this leads to an action 
of the renormalisation group on the space of defects. The fixed points of the renormalisation 
group are clearly of interest and these are the conformal defects. The problem of studying 
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conformal defects is equivalent to studying general conformal boundary conditions of a folded 
model [WA]. Even if one starts from a rational conformal field theory, the folded model 
will typically no longer be rational with respect to the diagonal symmetry preserved by 
the boundary condition corresponding to the defect. This means that the representation 
theoretic methods used to construct the bulk theory cannot be applied to get a handle on 
the defect itself. There are, however, two distinguished subsets of conformal defects which 
are known as topological (or purely transmitting) defects and factorising (or purely reflecting 
defects), which are much easier to study than the general case. 

Defects have an obvious generalisation to interfaces between different conformal field 
theories, possibly of different central charge. A particular model may, in fact, be simple 
enough to allow one to classify all conformal defects or interfaces; this is the case for defects 
in the Lee- Yang model |QKW 1 and in the critical Ising model [OA], and for interfaces between 
the Lee- Yang and the Ising model |]QRW |. Particular defects may preserve a rational sub- 
algebra in the folded model |QS| , |QRW| , |FG| , or a semi-classical analysis may suggest the 
existence of distinguished conformal defects and interfaces, for example in a WZW model 
for a given group [BG] or between such WZW models at different levels [FQ]. 

One may also think of interfaces as 'symmetries' which relate the properties of the two 
theories they link. For example, interfaces provide group symmetries and order-disorder 



dualities [Frl[| , they relate different renormalisation group flows of boundary conditions | GW ], 
there is a preferred interface joining the UV and IR fixed point of a given quantum field 
theory | BR2 ], and they can be used as spectrum generating symmetries in string theory 
[p3a2|1 . Interfaces can also be related to tunnelling in the quantum Hall effect [FFN]. Defects 
(and also interfaces) can be composed by placing the defect lines parallel to each other and 



letting their distance tend to zero [|PZ| , Fr2 , BRI , Ru, BB|. This process may or may not 
be singular, but if it can be defined then it describes an algebraic structure on the space of 
conformal field theories which deserves more investigation. 

In this paper we study a very simple example of a two-parameter space of renormalisation 
group flows. We identify the fixed points of these flows with the aim of finding new non- 
trivial conformal defects. The starting point for our flows is a particular topological defect 
D in the Virasoro minimal model M(p,p + 1) and the two parameters of our space are A; 
and A r , describing perturbations of D. We arrive at the following conjecture for the space 
of flows for p > 3 in the neighbourhood of this defect, shown in figure |]. 

Note that, in this figure, A; and A r represent renormalised coupling constants. An example 
of such coupling constants are the coupling constants used in the TCSA scheme in section [| 
taken at some fixed finite cut-off. 

Four of the fixed points (labelled D' and I) can be identified exactly and are topological 
defects; one can be studied numerically and we conjecture that it is a factorising defect F; 
finally we conjecture the existence of a new non-trivial conformal defect C which we study 
both numerically and perturbatively. 

A renormalisation group analysis in the case of WZW models showed a similar non-trivial 
fixed point for Xi = X r ^ 0, which is a good candidate for a non-top ological conformal defect 

S- 

The case p = 3 has recently been studied in QFFN ] and solved exactly by relating it to 
a free fermionic model. The space of RG flows is qualitatively different but the fixed points 
corresponding to /, C, D' and F are present in this model and agree with our conjectures for 
the exact forms of I, F and D'. Our numerical calculations agree with the results of [FFN], 
confirming the validity of our numerical method. 

Unfortunately, we have so far been unable to calculate exactly or numerically other 
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Figure 1: The proposed flows for the perturbation ( 2.23 ) for p > 3. The point D is the (1, 2)-defect. 
The possible endpoints are: I - the identity defect, D' - the (2, l)-defect, F - a factorising defect 
given by the sum X^r=i H r ' -0)(( r i 1 II °f P~ 1 conformal boundary conditions, and finally C - the new 
conformal defect. For details see the body of the paper. 



characteristic quantities of the new conformal defect C, such as its g- value |AL|| (defined 
to be that of the corresponding conformal boundary in the folded model) or its reflection 



coefficient [QRW]. However, we can calculate the g-value perturbatively for large values of 
p. 

This paper is organised as follows. In section ^ we collect the properties of topological 
defects needed in the subsequent analysis. Section |3] contains several exact results on defects, 
some new to this paper. Section [| contains the renormalisation group analysis of the pertur- 
bation for large p and identifies the new conformal defect as a perturbative fixed point. The 
numerial truncated conformal space approach used to support the proposed flows of figure 
H is described in section [5] with the results given in section |(| Finally, section ^ contains our 
conclusions. 



2 Topological, conformal and factorising defects 

To derive the properties of conformal defects, consider the complex plane with a defect line 
placed on the real axis. The defect preserves the conformal symmetry of the bulk theory if 
the field T xy = 4=(T — T) is continuous across the real axis. 

This condition is unchanged if we now consider the defect on a cylinder obtained by 
identifying z with z + 2tt. One can map this cylinder to the whole plane by z i— > exp(iz) 
and the real line gets mapped to the unit circle. In this formulation, a defect on the unit 
circle defines an operator D on the space of states H in radial quantisation. Translating 
the condition to this situation, a defect is conformal if the operator D commutes with the 
difference of the holomorphic and anti-holomorphic copy of the Virasoro modes, 

[L m - L_ m , D]=0 for all 77t G Z . (2.1) 

General conformal defects of a given conformal field theory are difficult to describe be- 
cause according to fl2,l|) they only preserve the diagonal Virasoro algebra = L m — L- m 
(with central charge c d = c + c) of the full holomorphic and anti-holomorphic symmetry. 
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Figure 2: The equivalence between a defect in a CFT and a boundary condition in the folded 
model 

A very useful way to think of this condition is in terms boundary conditions on the folded 
model. If one considers a CFT on the complex plane and folds the worldsheet over along the 
imaginary axis, the resulting model consists of the tensor product CFT <g> CFT on the half 
plane with a boundary condition inserted along the imaginary axis as shown in figure ^, where 
CFT is the original CFT with holomorphic and anti-holomorphic dependences swapped. If 



the original model had a conformal defect along this line, (2.1) implies that it is a conformally 
invariant boundary condition in the folded model. Note that this includes the possibility that 
there is no defect (or the identity defect) so that there is a distinguished conformal boundary 
condition corresponding to the absence of any defect. This correspondence between defects 
in a CFT and boundary conditions on the folded model CFT CFT will often be used in 
what follows. 

As said already, even if one starts from a rational model, the representation theoretic 
methods used to construct the bulk theory cannot be applied to get a handle on a general 
conformal defect. There are, however, two distinguished subsets of conformal defects which 
are much easier to study than the general case and which we consider now. 

2.1 Topological defects and defect fields 

Topological defects are a particular class of conformal defects which preserve a larger sym- 
metry than the diagonal Virasoro symmetry and are amenable to classification on account 
of this. The larger symmetry in question is, in fact, the full bulk symmetry, that is, the 
conformal defect obeys the stronger condition 

[L m , D] = = [L m , D] for all m G Z , (2.2) 

which clearly implies (|2.1| ). Such defects are called topological because they are tensionless 
and can be deformed on the surface without affecting the value of correlators. Topologi- 
cal defects were first studied in flPZfl (and the name 'topological' was introduced in [[BGfl ). 
Condition ( |2.2| ) is equivalent to demanding that T and T are separately continuous across 
the defect line. As a consequence, the space of defect fields forms a representation of the 
holomorphic and anti-holomorphic copy of the Virasoro algebra, just as the space of bulk 
fields. In particular, a defect field (ft has a left and right conformal weight (hi, h r ). 

The classification of topological defects can be carried out for rational conformal field 
theories and general modular invariants if one requires the defects to preserve the holomor- 



phic and anti-holomorphic copy of the chiral algebra [PZ, FRS |. Here we consider only the 
diagonal unitary Virasoro minimal models, for which the exposition simplifies. 

Let thus M(p,p+1) be the Virasoro minimal model of central charge c = 1 — G/(p 2 + p) 
and with diagonal modular invariant partition function. The irreducible representations Ri 
of the Virasoro algebra Vir occurring in these models are labelled by the set of Kac labels 



i G Ip, 



Tp = {(r, s) 1 < r < p— 1, 1 < s < p\/ ~ where (r, s) ~ (p—r,p+l—s) . (2.3) 
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The space of states propagating on a cylinder, or equivalently the space of bulk fields, 
decomposes into representations of Vir © Vir as 

H = Ri <g> Ri ■ (2.4) 



Both the elementary conformal boundary conditions [Ca] and the elementary topological 
defects [ |PZ| in M(p,p+1) are labelled by the set I p . Let ||o)) be the boundary state corre- 
sponding to removing the open unit disc from the complex plane and labelling the resulting 
boundary by a £l p , Denote by D k : TL — > TL the operator describing a topological defect 
with label k £ T p placed on the unit circle in the complex plane. Explicitly, these two 
quantities are given by 

x ""' and D k = ^ id Rilg)Ri , (2.5) 



where is the identity or vacuum representation (1,1), \i)) is the Ishibashi state in the 
(algebraic completion of) Ri (g) Ri and Sij is the modular 5-matrix for M(p,p+1), 

S (r , s) M = V« • • sin ■ sin ^ . (2.6) 

An important property of a defect is its entropy g defined in the same way as the boundary 
entropy and identical in value to the boundary entropy of the corresponding boundary con- 



dition in the folded model [A_L|. From ( |2.5| ) we read off the g- value of a boundary condition 
and a topological defect as 

9(\\a))) = ^$= and g(D k ) = f° . (2.7) 

We recall that topological defects can be deformed freely on the surface as long as they 
do not cross field insertions, boundaries, or defect lines. The fusion of two topological defects 
corresponds to the composition of the defect operators, and the fusion of a defect with a 
boundary condition is given by the action of the defect operator on the boundary state. One 
easily checks that for i,j, fc,o£ I p , 

Z> (lil) = id n , D k D l= J2 N kT An , D k || a)) = Y, N ka m , (2-8) 

where are the fusion rule coefficients. By evaluating the corresponding partition func- 
tions, or by using the methods of |[Fr2| , one finds that the space of defect fields living on a 
defect labelled by k G X v is given by 

H%= (Ri^Rj)®^ 1 ?^^ . (2.9) 
i,jei P 

We will also need the space of states propagating on a strip with boundary condition (1,1) 
on one side and a £ I p on the other side, and with a defect labelled k &X P running parallel 
to the boundaries inside the strip. In the same way as ( |2.9| ) one finds this is 

H (W = 0^ . (2 .io) 
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Figure 3: The equivalence between a factorising defect in a CFT and separate boundary 
conditions in the folded model 



2.2 Factorising defects 

By a factorising defect we mean a conformal defect represented by an operator F that satisfies 
the stronger conditions 

(L m - L_ m )F = = F{L m - Z_ m ) for all m € Z , (2.11) 



which then imply ([2.1|). Such a defect is totally reflecting and is simply a sum of products 
of conformal boundary states 

F= n ^H a »« 6 H ( 2 - 12 ) 

for some non-negative integers n ab . In the folded model this corresponds to separate bound- 
ary conditions on the two sheets, as in figure ||. The g- value of the factorising defect ( 2.12 ) 
is 

9(F) = £ n ab ^ . (2.13) 



2.3 Conformal defects 

Returning to general conformal defects, we give the form of the partition function on a 
cylinder and define a useful pairing on such defects. 

Let d and d' label two conformal defects and consider a torus of width R and height L 
with two vertical defect lines at distance r/R, < rj < 1, from each other. Denote by D d and 
Dji the defect operators of the defects d and d'. Define the partition function 

Z d>d , (q; r,) = tv n ( qViLo+lo-c/12) ^ ~(l- v )(L 0+ l -c/l2) ^ ^ (2 M) 

where q = e~ 2nL / R and q = e ~ 2nR / L . The adjoint (D^ describes the defect operator for the 
defect d inserted with reversed orientation (see [ |PZ| ] and |Fu| , sect. 6.1]). Using the explicit 
expressions ( |2.5| ) one checks that for the factorising and topological defects in M(p,p+1) one 
has 

(||o))((6||)+ = ||6))((a|| and = D k . (2.15) 

In the special case that 77 = ^, we can identify the original CFT on a torus with two 
defects with the folded model on a cylinder of circumference L and width R/2 with conformal 
boundary conditions corresponding to the defects d and d' at the ends of the cylinder. The 
partition function of the folded model on a cylinder of circumference L and width R/2 is 
a sum of Virasoro characters Xh,2c{<i) at twice the central charge of M(p,p+1), so that the 
partition function can be written as 

Zd,d'{m 2) = ^rnhXhiciq) (2.16) 
h 

where q = e ~' 2nL / R i mh \ s the multiplicity of the representation with highest weight h and 
the sum can be over a finite or infinite set of weights h. We have assumed that the spectrum 
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of the folded model on the strip is discrete, i.e. that Z^d'{.<l]\) is a sum of powers of q, 
rather than an integral. This may not be the case for a general conformal defect, but it is 
true for the topological defects and for the totally factorising defects ||a))((6|| as can be 



checked with the explicit expression in (2.5). We investigate defects which can be reached 
by renormalisation group flows starting at these defects, and so we expect the spectrum at 
the IR fixed point to be discrete as well. 
We define the pairing (d, d!) by 

(d,d') =m G Z> , (2.17) 

i.e. the multiplicity of the vacuum character in Z^ ^^q; ^). For unitary models conformal 
weights have to be non-negative so that we can write (d, d') as the limit 

(d,df) = limcf/ 12 Z d4 ,{q-\) . (2.18) 

Since tr(X) = (tr(Xt))* anc i (d,d') is real, it follows from ( fTPf ) that (d,d') = (d',d). 

By passing to the folded model one checks that (d, d) > 1, that (d,d) = 1 if and only 
if the conformal defect d is elementary, and that (d, d') = if d and d' are elementary and 
distinct. The latter statement follows in the folded model since one cannot have a boundary 
changing field of weight zero between two distinct elementary boundary conditions a and b 
(such a field could be pushed along the a boundary changing it into a b boundary without 
affecting correlators, showing that a and b have equal boundary states^ ). 

In particular, if d is a superposition and d' is elementary, then (d, d') gives the multiplicity 
of d' in the decomposition of d into elementary defects. 

Another useful property of the pairing is that if k € I p is a topological defect, and k * d 
denotes the fused defect (which is again conformal and has the defect operator D^Dd) we 
have 

(k*d,d') = (d,k*d') and (d ★ k, d!) = (d, d' * k) . (2.19) 
The first equality follows from 

tr W ( q^o+Lo-c/U) ( DkDd) \ ^(Lo+Lo-c/12) ^ 

(2.20) 

tr w ( qi(L +L ~c/l2) {Dd)] ~A(Lo+Lo-c/i2) DkDd/ ■ 



as -Dfc is self-adjoint and commutes with Lq + Lq. The second equality can be seen similarly. 
2.4 Defect perturbations 

Finally we consider perturbations of topological defects. A perturbation of a topological de- 
fect by a defect field (t>h u h r is relevant if h[ + h r < 1. This is in contrast to bulk perturbations, 
which are integrated over the whole surface, not just over the defect line, and are relevant if 
hi + h r < 2. 

The first thing to note is that all topological defects in unitary minimal models have 
relevant perturbations, even the identity defect (that is, no defect), which is quite different 



Some care has to be taken since in general the boundary state does not characterise the boundary 
condition uniquely. Similarly, the defect operator does in general not determine the defect uniquely. An 



example of this is discussed in [Fu . For unitary models this ambiguity is not expected to occur, i.e. there is 



a one-to-one correspondence between boundary conditions and boundary states, and similar for defects. 
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from the boundary situation where there are (p— 1) stable boundary conditions for M(p,p+1). 
The fields on the identity defect are all scalars and are exactly the bulk fields 4>h,h- There 
are \_\j2p 2 + 2p + lj — 2 of these with scale dimension less than or equal to 1, as opposed to 
2p — 3 which have scale dimension < 2. In any case, even the identity defect is unstable to 
defect perturbations by these bulk fields (considered as fields on the defect). 

As another example, on the (1,2) defect there are both scalar primary fields §h,h an d 
non-scalar primary fields 4>h,h' with h h! . There are (p — l) 2 scalar primary fields of which 
2[y/2p 2 + 2p + lj — 6 have dimension < 1 and (p — l)(p — 2) non-scalar primary fields of 
which 2{p — 2) have dimension < 1. The relevant non-scalar fields are 

{(f>(r,r),(r,r+2) and <t>(r,r+2),(r,r) | 1 < r < p — 2} . (2.21) 

We note here that the interpretation of a defect as a boundary condition in a folded 



model means that the boundary ^-theorem [AL, FK] applies to perturbed defects as well: g 



is decreasing along RG flows and so restricts the possible IR fixed points accessible from any 
given UV fixed point. 

The simplest class of relevant perturbations are those by chiral defect fields (f>ho ° r 4>o,h 



with h < 1. Such perturbations were investigated in [Ru]. These perturbed defects are 



believed to have particularly nice properties, for example to commute with Lq + Lq, see 



section 2.4.1 below. 

On the other hand, since a topological defect perturbed by A; • <ph,o, say, will commute 
with all L m for arbitrary values of the coupling A;, the conformal defect obtained as the IR 
fixed point for large A; will necessarily be topological [JBGJ]. To obtain an IR fixed point 
which is conformal but not topological, the next simplest perturbation to try is the defect 
field 

A* • </> h , + A r • (j) ,h ■ (2.22) 

A topological defect perturbed by this field will no longer commute with either the -L m 's or 
the L^s and in general it will also not commute with Lq + Lq (there are exceptions, for 
example the fusion of a defect perturbed by (ph,o with a defect perturbed by (/>o,/i). 



For the conformal weight h in ( 2.22 ) we choose ^(1,3) = (p — l)/(p + 1). This field is 
of interest for several related reasons: in many circumstances it leads to integrable pertur- 
bations which can be solved exactly and we expect the purely chiral perturbations to be 
integrable and exactly solvable; this perturbation is also related to one of the integrable 
lattice description of minimal models; finally, it is particularly suitable for a renormalisation 
group analysis, as it does not generate further relevant fields under fusion (note that <ph^ 
already has weight greater than or equal to 1), and it becomes marginal in the limit p — > 00. 

The simplest topological defect which allows for a perturbation by ( |2.22| ) is the (1,2)- 
defect, and so this is finally the situation which we will study: 



D l2 {\i(f> + \ r <j>) 

The (l,2)-defect perturbed by Xi<p + X r <j) with 4> = 4>h a 3) ,o and 



(1,3) i u V VU,/l( lj3 )- 



(2.23) 



There are two useful results which we can derive for this perturbed defect which we give 
now. 

2.4.1 Translation invariance 

As mentioned above, it is believed that a topological defect perturbed by a chiral field still 
commutes with Lq + Lq. This is shown in [|Ru|| for chiral fields with h < 1/2 where there 
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are no UV divergences in the expansion of the perturbed defect; for fields with h > 1/2 this 
property depends on the existence of a suitable regulator which preserves the commutators. 
Assuming for the moment that such a regulator can be found, in the present case this would 
imply 

[L o + Z o ,£>(i,2)(A0)] =0 and [L + L Q , D (1 , 2 )(A0)] = . (2.24) 

This guarantees that the action of -D( 12 )(A</>) and -D( 12 )(A</>) on a boundary state is non- 
singular. 

One particular consequence of the property (|2.24j) is that the spectrum of states of the 
model on a strip containing a defect parallel to the edges of the strip is independent of the 
position of the defect. This is easy to see in the case of the TCSA regulator, as we prove in 



appendix A. 4, so that it is very reasonable to apply this assumption to our TCSA results 



and we will use this in section [3T] to deduce the end-points of the chiral perturbations. 
2.4.2 Commutation with a subset of topological defects 

Let D(Yl a A«0«) be a topological defect perturbed by defect fields 4> a where cj) a is a field 
(not necessarily primary) in the sector Rk a (g> R[ a . Suppose that either all of the k a are of 
the form (1, s a ) with s a odd, or all of the l a are of the form (1, s a ) with s a odd. Then 

[%),O(E tt U)]=0 farl<r<p-l. (2.25) 

We will demonstrate this in the case that the k a are of the form (1, s a ) with s a odd. The 
second case can be seen analogously. 

Let x € Ri <S> Ri and y £ Rj <S> Rj for i = (a, b) and j = (c, d) elements of X p . From the 



explicit form of Dr r> -n in (2.5) and abbreviating D = D(^2 a \ a <fia) we find 



(x\D {rA) D\y)= S{r ' 1){a ' b) (x\D\y) , (x\D D^y) = ^ {x \ D \ y) , (2 .26) 
a (l,l)(a,b) a (l,l)(c,d) 

Now expand out the exponential in (x\D\y). Since each term in the integrand is in the sector 
Rn iSa ) <8> Ri a , the expression (x\D\y) can be non-zero only if i and j are such that a=c and 
b—d is even. Substituting the expression (|2.6|) for S'y gives 

S (r,i)(a,b) = ( -_ 1 ^(r+l)b sin(7r(p+l)ra/p) 
s (i,i)(a,b) sin(vr(p+l)a/p)) 

Since a=c and b—d is even, the two expressions in (|2.26| ) are therefore equal for all states 
x,y, proving ( |2.25| ). 

For us the commutator (|2.25|) is interesting because it implies [ D( r ^ , D( 1 2)(Xi(f> + A r ^)] = 
with (j) and 4> as in ( p. 23 ). However, let us note in passing that (2.25) can also be applied 



to bulk perturbations. In particular, we can set D = exp(— £H) with H the Hamiltonian of 
the minimal model perturbed by the bulk field <j>(\$\(\$\ since 

H = J ( U T + f ) + ^(l,3),(l,3) ) dx (2.28) 

and T and T are descendent fields in the (1, 1) representation. The topological defects -D( r ,i) 
continue to commute with D for all values of the coupling /x, and thus they are conserved 
charges also off criticality. This provides an independent argument for the observation made 



in | FGS | that under the bulk flow between two minimal models generated by <t>(x^)Ax^)i the 



topological defect -D( r ,i) should again flow to a topological defect. 
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3 Exact results 



The properties (|2.24|) and (|2.25| ) of -D( li2 )(A/0 + K4>) enable us to find two exact results 



for the perturbed defect, which we describe in sections |3.1| and |3.2| . We also summarise in 
section |3.3| the known results for the Ising model (the case p = 3) which is distinct from the 
other models. 

3.1 The chiral perturbations 

The translation invariance of the purely chiral and purely anti-chiral perturbations in (2.24) 
allows us to determine the nature of these flows exactly. 



Consider the situation described in ( 2.23j ), i.e. -D(i,2) (\i<p + \4>) is the operator obtained 



by placing the (l,2)-defect on the unit circle and perturbing by A/0 + A r </> with A; and 
A r G IR. Let t/j be the primary boundary field on the (1, 2)-boundary condition in the -R(i,3) 
representation and denote by || (1, 2) + Xip}} the boundary state obtained by perturbing the 
(1, 2)-boundary by Xip. One can normalise the fields (f>, 4> and ip such that the identity^ 

D(i,2)(m\ 1, 1)) = II (1, 2) + AVO) = £>(i,2)(A0)|| 1, 1)) (3.1) 

holds. The boundary flows in the middle of this series of equalities are already known 
[ LSSj [RRSj , |GRWf , 



1,1)) , -°°- A |] (1,2) H- AV» X ~* +C ° > ||2,1». (3.2) 



As mentioned in section 2.4 , a chirally perturbed topological defect necessarily has again a 
topological defect as an IR fixed point. Acting with the IR topological defect on || 1, 1)) gives 
a conformal boundary state that has to agree with ([T^). The explicit expressions in ( |2.5[) 
show that a topological defect is uniquely determined by its action on ||1, 1)), so that the 
above reasoning fixes the horizontal and vertical flows in figure |l| to be 



/ = ■ %, 2 ) (k<f>) ■ D' = -0(2,1) 

I = D( lt y < r D( lt2 )(\ r <t>) — r ~ >+ °° > -D' = £>(2,i) 

3.2 The factorising component of any IR fixed point 



(3.3) 



The fact that our perturbed defects commute with the set of topological defects {D^ t ^} 
shown in ( |2.25| ) allows us to determine the building blocks of any factorising component of 
any IR fixed point. 

Consider the superposition of elementary factorising defects given by 

„<,!»' „ „ ;,, „ \h ■ V odd and s = s' = 2±i 

^=E^Dii.'»«iyi^ q -{; : l erwise • 

The factor of \ has to be included because for s = s' = ^i- we have \\t, ^-)){{t, 2i- || = 
\\p — t, ^t-))((p — t, ||, so that in this case each factorising defect appears twice in the 
sum. For example, in the Ising model (p = 3) we have F 2 ' 2 = || 1, 2)) (( 1, 2 1| . 



2 This identity is evident in the case A = 0. Since the space of (1, 3)-primary fields on the (1, 2)-boundary 
is one-dimensional and spanned by ip, the perturbing fields <j> and 4> have to be proportional to ip once the 
(l,2)-defect is p ushed on top of the (1, l)-boundary. Alternatively, the identity can be proved with the 
methods in lPr2 
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It is shown in appendix A.l that [Dr ri i\, F s \ s '] = for r = l,...,p— 1, and that the 



following statement is true: If a conformal defect C obeys [Dm i), C] = (which implies that 
[Dr r u, C] = for all r) then we can write 

C = R + F (3.5) 

where R is a conformal defect that does not contain factorising defects as summands, i.e. 
(R, \\a)){{b\\) = for all a, b G Z p , and where F is a combination of factorising defects taking 
the form 

F = E E<' S,jD (M)^ |S ' (3-6) 

s,s'=l a=l 

for suitable constants G Z>o- If C is an IR fixed point of D(i 2){\i4> + ^r4>) in a 

particular direction in the (A;, A r )-plane, it will commute with the -D( r ,i) and hence can be 
written in the form ( p.5| ). This will be exploited in the (^-function analysis in section |4| and 
in identifying the spectrum on a strip with a perturbed defect via TCSA in section |(| 

3.3 The Ising case 

The Ising model is unique amongst unitary minimal models in that its conformal defects can 
be completely classified. The reason is that its square Ising<g>Ising is a c = 1 CFT which can 
be described in terms of an orbifolded free boson. Using this relation, Oshikawa and Affleck 



give a complete classification in [ OA ] of the conformal defects of the Ising model in terms of 
the conformal boundary conditions of the orbifolded boson as follows. 

There is a continuous family of Dirichlet boundary conditions D(tpo) subject to -D(yo) = 
D(—ipo) = D((po+27r) for which a fundamental domain is c^o £ [0, 7r]; and a continuous family 
of Neumann boundary conditions N(ipo) with N(<po) = N(—(Pq) = N((f>Q + ir) for which a 
fundamental domain is <p G [0, 5]. The boundary conditions in the two continuous families 
are all elementary, except for the boundary conditions at the endpoints of the fundamental 
domains, which each split into two elementary boundary conditions so giving a discrete set 
of eight boundary conditions. Altogether we have 

\D((f )) with <p G (0,tt) and \D(0)) ± , \D(ir)) ± , 
\N{<p )) with O G (0,tt/2) and \N(0)) ± , |A^(vr/2)) ± . 

These have g values y/2 for the continuous Neumann series, 1 for the continuous Dirichlet 
series, l/v2 for the discrete Neumann defects and 1/2 for the discrete Dirichlet defects. 

There are three topological defects in the Ising model given by Kac labels (1, 1) = id, 
(1, 2) = a and (2, 1) = e and three Cardy boundary conditions "+"= (1, 1), "/"= (1, 2) and 
"— "= (2, 1). The topological and factorising defects are identified with the orbifolded boson 
boundary conditions [OA, QRW| ] as shown in figure ||. 



To identify the fixed points of the numerical RG flows in section |6.l| we need the partition 
function encoding the spectrum of states on a strip with the (1, 1) ="+" boundary condition 
on the two edges and a defect in the middle. When this strip is folded, this is equal to the 
partition function in the orbifolded boson model with the (++) = -0(0)+ boundary condition 
at one side of the folded strip and some boundary condition corresponding to the Ising defect 
at the other boundary. In the UV this is the (1,2) = a = N(ir/4:) boundary condition, but 
by the (/-theorem, the IR point fixed point can be any Dirichlet boundary condition (together 



12 



++ 



D(0) + . 1 = ^(1) e = -D(^) •(-+) = £>(tt) + 



(— ) = D(0)_ • (//) = • (+-) = D(tt)_ 

(/+) = AT(0)+ • <r = -|V(|) • (+/) = iV(f )+ 

(/-) = AT(0)_ • ' • (-/) = N{1 )- 

Figure 4: The Dirichlet (top) and Neumann (bottom) boundary conditions on the orbifolded 
free boson and the corresponding topological and factorising defects in the folded Ising model 

with any of the four discrete Neumann conditions, but these do not arise in the flows we 



consider); the partition function of such a system is given in [OA] 



' n=—oo 

where q = e - 2nL / R . This expression is also valid at the endpoints (fo = 0, tt if one takes 
D(0) = D(0) + + D(0)_ and D(ir) = D(ir) + + D(n)„. 

As has been discussed in section |3.1|, the IR fixed points of the purely chiral flows can 
be easily found: they are the (2, 1) = e and (1, 1) = id topological defects in the positive 
and negative directions, respectively. The remaining fixed points have been recently found 
(amongst other results) in [FFN] and the space of flows is shown in figure |l~3|a. In that paper, 



the perturbed defects are parametrised by an angle 6 with the general result 

AT(^o) — > D(<p + 6) . (3.8) 

We define the angle a by tan(a) = A r /A;, with a G (— f , f ) for A; > and a G (^, ^) for 
A; < 0. One can check that a is related to 9 by a = tt — 9 so that the prediction of [ ]FFN ] 
for the RG flows of the defects ( [2.23D is 

X (7 = N{tt/A) — ► £>(57r/4-a), (3.9) 

which is exactly what we shall find in section |6.1|. 



4 Perturbative analysis 

The renormalisation group equations for perturbations of conformal boundary conditions 
(and hence also of defects) have been known for a long time, studied first in [|AL| and used 
in |RRS , Gr] to study the flows in unitary minimal models. 



The necessary ingredients are a conformal boundary condition, a set of relevant boundary 
fields S = {4>i} which is closed in the sense that all the relevant fields occurring in the fusion 
of any two fields in S is also in S, the conformal dimensions hi of these fields, and the 
structure constants o- L jk between these relevant fields. Furthermore, since the perturbative 
integrals are divergent for hi > 1/2, it is necessary to regularise them to which end we 



introduce a UV cut-off a as in [AL] and define dimensionless couplings /Uj = \ia y% where 
yi = 1 - hi. 
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Using the regularisation in [AL], the RG equations for the couplings fi{ are 



fj-i = Ui^i ~ ^2 CijkfJ-j^k + 0(fi 3 



(4.1) 



The change in the g- value of the boundary condition/defect is also known in perturbation 



theory | AL , R.RS ] and is (to third order) 
log(ff(Mi)) - log (go) = 



(4.2) 



The scale dimensions of the perturbing fields at a fixed point can be read off by linearising 
the RG equations about the new fixed point, the eigenvalues being 1—h' where h' is the scale 



dimension of the field [ Gr ] . 

We now restrict attention to the perturbed defects ( 2.23| ). In this case the space of 
relevant defect fields generated by <f> and 4> closes on those fields alone, that is they generate 
no new relevant defect fields, so we do not need to introduce couplings to any other fields 
into the RG equations. 

Secondly, since <j> and <j) are chiral and anti-chiral respectively, the only structure constants 
appearing in the RG equations which can be non-zerof] are and c^^, and these are in 
fact equal and the same as the boundary structure constant cy,^, for the field ip = ^(1,3) on 
the (1, 2) conformal boundary condition (see footnote [2] and [Ru, sect. 2]). This means that 
the defect RG equations are two decoupled copies of the boundary RG equations of [RRS], 
that is 



fir 



(4.3) 



'1,3 = ^pr ana - C V>V>V> 
given to third order by 



where y = 1 — h\ 



8/3 + 0(y). The g value of the perturbed defect is 



log(g(fli, fJL r )) - log(flr(0,0)) 



-TT 2 y(fif + fi 2 r ] 



27r2 / 3 3' 



(4.4) 



fi* (where fi* = y/c^^) so that the 



Each RG equation has two fixed points, fi = and fi 
combined RG equations have four fixed points. 

It is important to note that the perturbation expansion is for small y, that is for large 
p where the central charge is close to 1. While the perturbative results for boundary flows 

3, this will not be the case here, as the 
3. 



in [RRS] still hold true all the way down to p 
spectrum of the IR fixed points is quite different at p 
The perturbative RG flows are shown in figure ||| 

The values of the couplings at the fixed points, the corresponding value of g and the scale 
dimensions of the perturbing fields at the fixed point (to first order in y) are given below 



(W,Mr) 


log {9 /go) 




h$ 


(0,0) 





1-2/ 


1-2/ 


(/AO) 


ttV 

2 8 3 


1 + 2/ 


1-2/ 


(0,/i*) 


Iff 
2 8 3 


1-2/ 


1 + 2/ 


(/Am*) 


rf 
4 


1 + 2/ 


1 + 2/ 



(4.5) 



It is important to note that c^^,^ 7^ for p > 3 but when p = 3 the constant c^^^ vanishes, so that this 
case is not directly amenable to the analysis here. 
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Figure 5: The perturbative flows for the system ( |4.3] ). The perturbative fixed points are as 
in figure |], the black diamond being the conformal defect C. 



The topological defects of lowest g value and their expansions in y to third order are 



Defect 


9 


log(g) 


^(i,i)' d {p-i,i) 


1 





D (2,l)i D (p-2,1) 


2cos(£) 


log(2) - _ 


D (l,2)i D (l,p-l) 




log(2) - ^ 


D (3,l), D (p-3,l) 


l + 2cos(f ) 


9 9 9 ^ 

log(3) _ ZE^L _ 2lM 




l + 2cos(^ r ) 


log(3) - ^ 



(4.6) 



Comparing 5- values and noting that the end point of an RG flow generated by a chiral field 
must be topological, we see that the chiral and anti-chiral perturbations of -D(i,2) can a priori 
have either the -D(2,i) or ^(p-2,1) defect as their perturbative fixed point. This can be decided 
by considering two point functions of bulk fields in the presence of the defect, which can be 
calculated using fl2.5|) : 





D(l,2) 


-0(2,1) -D(p-2,1) 


(¥'(12,12) 
(^(21,21) 


D\<f(l2,12)) 
D \<P(21,2l)) 


-2cos^ 
- 2cos ^T 


-2cosf -2(-l)fcos| 
-2cos^ -2(-l)f +1 cos^ 



(4.7) 



We see that the two-point functions in the presence of -0(2, 1) are perturbatively close to 
those in the presence of -D(i ; 2)) whereas in the presence of -D( p -2,i) the two point functions 
of y?(i2,i2) anci ^(21,21) are n °t perturbatively close for p odd and p even respectively. Thus 
we deduce that the perturbative fixed point for the chiral perturbations of Dn 2 ) are -D(2,i) 
and not -D( p _2,i)- 

Equally, we see that the C defect which is the fourth fixed point in the '++' direction is 
not a topological defect. The question remains whether C can be found as a superposition 
of factorising defects or the Identity defect superposed with a set of factorising defects. 

This is easy to address given the general form of a fixed point in (|3.5D . We first compute 
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the g- value of each summand in (|3.6| ) 



a(D ^ (s^fl) ( S ^M\ ( S ^A ( n^^M) 5 

9{ (0,1) J " ) ' ™(f) J 

= a^'^-^^+O^ 2 ). (4.8) 

Here 5 is either | or 1 as in ( |3.4j) ; it does not appear in the second line as this is evaluated 
for fixed s, s' and large p. 

Since each of the combinations D( a ^F s \ s ' has a g- value greater than or equal to one 
as does each of the topological defects, and since the g-value of -0(1,2) is l ess than two, 
it is clear that any IR fixed point of a flow starting from D(i,2) can include at most one 
topological defect or one factorising combination but not both. Furthermore, for large p, the 
only possible factorising combinations with small enough g value are F 1 ! 1 and D^^F 1 ^ 
but these do not agree with the perturbative calculation of the g-value of C. 

(We note here that there is an additional candidate in the case p = 4: the factorising 
defect D(2 a ^ so nas a S-value lower than that of Dn 2 ))- 

Summarising, from an analysis of (/-values, C is not a linear superposition of topological 
defects and factorising defects and so contains a new conformal defect; the most likely result 
is that C is a new elementary conformal defect, but that cannot be determined from the 
calculations of the g value alone. 

There remains the question of identifying the nature of the fields eft and 4> at the fixed 
points. We can do this by checking which fields on D(2,i) have weights as given in table ( ^4.5[ ). 
This is straightforward and the only candidates are as given below 



(A»i>A*r) 


D 4> (j) 


(0,0) 


£>(1,2) 


0(13)(11) 


0(11)(13) 


(/AO) 


-0(2,1) 


0(31)(11) 


< / > (33)(13) 




£(2,1) 


( / > (13)(33) 


0(11)(31) 




c 


? 


? 



This means that we should be able to obtain the 'C defect by TCSA calculations starting 
both from -D(i,2) an d perturbing by <f> + 4> an d also starting from £(2,1) an d perturbing by 
0(13,33) ( or equivalently by 0(33,13))- This is one of the subjects of section ||. 



5 Truncated conformal space approach for defects 
5.1 The TCSA Hamiltonian for defects 

The so-called truncated conformal space approach or TCSA is a numerical method for cal- 
culating the spectrum of a perturbed conformal field theory. One chooses a system which 
admits a Hamiltonian description and then the Hamiltonian is restricted to a finite dimen- 
sional subspace of the infinite dimensional Hilbert space spanned by energy eigenstates whose 
energy eigenvalues (or conformal weights) are not greater than a threshold value. The orig- 



inal idea was proposed in YZ] and it was applied for the first time for boundary problems 



in I Do 



One can envisage several possible systems involving defects which allow such a Hamil- 



tonian description. We choose to apply a generalisation of the boundary TCSA of [Do], 
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Figure 6: The strip is mapped to the upper half- plane and a defect running parallel to the edges of 
the strip is mapped into a defect running in a radial straight line. The dotted line shows an equal 
time slice for the Hamiltonian description. 



because the space of states is much smaller than in the other possibilities (they form irre- 
ducible representations of a single copy of the Virasoro algebra - see (|2.10D ). This makes 
the numerical calculations easier and we can also make use of the known properties of the 
interactions between defects and boundaries. The simplest situation is given by a strip of 
width R with conformal boundary condition (1,1) on both sides and a defect labelled by 
k G Ip running parallel to the boundaries at a distance a from one boundary. According to 
( |2.10| ) the Hilbert space consists of a sole representation (independently of a), 



n 



(5.1) 



By the exponential map the strip unfolds onto the upper half plane with the boundaries 
lying along the negative and the positive real axis and the defect line running from the origin 
to infinity at an angle 6 = an/R, as shown in figure ^. 

The Hamiltonian of this conformal field theory perturbed by ./V defect fields is 



H 



7T 

R 



L, 



c 

24 



JV 



i=i 



IT 



(5.2) 



We use a non-orthonormal basis consisting of vectors of the form 



L-m ■ ■ ■ L- 



m > ■ ■ ■ > n m > . 



(5.3) 



We label these basis vectors \vi), their Lq eigenvalues by Aj and their inner product matrix 
Gij = (vi\vj). The inner product matrix appears explicitly in the matrix to be diagonalised: 



' - 1 



N 



i{h\-K) 



(5.4) 



where we have introduced an energy scale A so that the operator h 
and its matrix elements are hij where h\vj 



H/A is dimensionless 



hij\vi) 



furthermore, B{0) l - k are the matrix 
elements (vj\4> h i hi (e t0 , e~ l0 )\vk) in the basis (|5.3|), r = AR is a dimensionless parameter and 

Ki = Xi(R/7r) 1 ~ h ^~ h ' r are dimensionless coupling constants. 

In the general case it is necessary to determine the structure constants to normalise the 
matrix elements B{9) % correctly, but in the cases we consider the Hilbert space consists of 
a single representation so that there is a single highest weight \hk) and so we normalise the 
fields so that {h k \<f> h .^ ,(l)\h k ) = 1. 
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The two special cases we will consider are the case N = 1 of a single non-chiral pertur- 



bation and the case N = 2 where the perturbing fields are two chiral perturbations (2.22) 
with (hj,hl) = (h,0) and (hf,h%) = (0,h). In the first case, there is a single dimensionless 
coupling constant k and the Hamiltonian becomes 

r - h jk = (Aj - ^)6 jk + k jfo-h,)9 {G -i B (e)) jk . (5.5) 

In the second case, the two coupling constants are K\ = K\ = kcos(o) and K r = k>2 = Ksin(a) 
so that 



h jk = (A, - ±)8 jk + k (cos(a) e* hd {G~ l B{9)) jk + sin(a) e~ ihe {G~ 1 B{9)*)j k ) , (5.6) 



where B(9)* is the complex conjugate of the matrix B{6). 
From the equality 

(A;- A^OH^OsOh) = («i|[£o 5 0h(*)]|«j) = ( h + z §^j H<K*)h> . ( 5 - 7 ) 

it follows that the matrix elements of a chiral field between Lq eigenstates have coordinate 
dependence 

(v i \ct> h {z)\v j ) = C-z^- h -^. (5.8) 

This implies that 

e ihe B(8) M + e- ihe B(d)* kl oc e*^*"^^ + e^ 1 '^ 6 . (5.9) 

When the defect is in the middle of the strip, 9 = tt/2, and the couplings to the left and 
right fields are the same (cos(a) = sin(a)) then the (k,l) matrix element of the perturbing 
Hamiltonian is proportional to i^-^i _|_ ^_^A(-A fc _ j£ ^ _ ^ - g an 0( j^ integer, the two 

terms cancel each other and the matrix element becomes exactly zero. In all our examples, 
the Hilbert space consists of a single representation, so the weight differences are always 
integers and they are odd if one of the vectors belongs to an even level while the other one 
belongs to an odd level. This means that for 9 = tt/2 and Ki = K r , the even and odd levels 
of the representation constitute two disconnected sectors: the Hamiltonian has zero matrix 
elements between vectors belonging to different sectors and thus it is block diagonal. 

5.2 Finite-Size scaling in TCSA 

There are two RG-type flows at play in the TCSA system. The first is the physical flow 
in which we are interested, that is the flow induced by scaling the size of the system. We 
parametrise this flow by a parameter t. The second is that induced by changing the cutoff 
iV limiting the size of the space of states. We are interested in the case of two relevant fields 
with a Hamiltonian (|5.6| ) in which case there are three couplings to consider, r(N,t) and 

Ki(N,t). 

Considering the case of infinite cut-off first, we expect that the TCSA system approaches 
the system with canonical scaling, so that at N = oo 

r = r r = rne* 

yt ( 5 - 10 ) 

Ki = yKi Ki = k\ e yi 



18 



For finite N, these are both altered. This has two effects: firstly, the eigenvalues of hjk are 
related to those of the renormalised Hamiltonian up to an unknown factor, so that it is only 
ratios of energy differences that can be calculated accurately; secondly, the TCSA flow can 
approximate IR fixed points for finite values of the bare coupling constants or equivalently a 
finite value of the bare volume [|Fe| , so that we can find one, or indeed more, conformal field 
points for increasing values of k. 



We can find the change in the coupling constants as N changes, which from Fe] is 
according to an RG-type equation of the form 

-N^(N,t) = h{<N,t) ] N) , (5.11) 

where 

(kfaN) = -c ijkKj K k N~ y + o{k 2 ,N-\ V ) . (5.12) 

We can use this equation to find the finite-size scaling flow at fixed cut-off under two as- 
sumptions. 

The first assumption, already stated, is that in the N —* oo limit, the TCSA couplings 
become the "linear beta function" couplings, a fact assumed in much of the literature, so 
that the finite-size flow at N = oo is 

«*(«>,*) = e yt . (5.13) 

The second assumption is that the leading terms in the beta-functions take a simple scaling 
form, that is 

AO = Ny 7i (KN-y) + . . . . (5.14) 

This is true for the quadratic terms and a reasonable assumption for the leading behaviour 
of the higher order terms. 

Given these two assumptions, we can then deduce the following three results (see ap- 



pendix A. 2 ) 



Ki{Ne\t) = e»'/i(« ) (5.15) 
Ki(N,t) = e yt Ki(JVe"*,0) (5.16) 

^(N,t) = f3i(K(N,t);N) (5.17) 

The first result states that we can follow a finite-size scaling flow by scaling the couplings 
by the same factor but also increasing the cut-off. This is useful if we do not know the 
beta-functions, as is the case here. The second result gives the scaling flow at fixed N in 
terms of the couplings at smaller N, and the third result states that the finite-size flow at a 
fixed cut-off is governed by a standard beta-function relation, where the functions are simply 
related to those for the change in N: 

PifaN) =yKi+Pi{K;N) . (5.18) 

Finally, a conformal defect corresponds to a zero of the beta- functions @i(k;N). From the 
leading behaviour ( 5.14j ), the positions of the zeroes n*a(N) scale approximately with N as 



N v which we check in one case in section y. 

Given the TCSA Hamiltonian, it is possible to calculate the TCSA approximation of 
many different physical quantities, but in this paper we have examined only one, namely the 
spectrum of the Hamiltonian, and we describe in the next section how we use this to identify 
the presence of possible RG fixed points, that is conformal defects, and attempt to identify 
them. 
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5.3 Identification of conformal defects using TCSA 

From the TCSA method, we obtain the spectrum of the Hamiltonian on a strip. The simplest 
way to identify candidate conformal defects is to place the defect in the middle of the strip 
and use the fact that in this case the partition function takes the special form ( p. 16 ) and 



the spectrum is organised into representations of the Virasoro algebra of central charge 2c. 
Consequently we examine the TCSA results to find spectra which take this special form and 
propose these as candidates for conformal defects. 

To be explicit, we consider the model on a strip of width R with conformal boundary 
conditions on the two edges and the defect D in the middle. In our TCSA calculations we 
will always take the boundary condition to be the same on both edges, but for the moment 
we shall label them B and B' for generality. This system can be folded into the folded model 
of central charge 2c on a strip of width R/2 with a factorised boundary condition (B, B') on 
one edge and a boundary condition corresponding to the defect on the other edge. If this 
defect is conformal, then the spectrum falls into representations of the Virasoro algebra of 
charge 2c and 

Z (B,B');D conformal =^2rn h Xh,2c(l) > ( 5 - 19 ) 
h 

where q = exp(— 2ttL/R). This means that the spectrum falls into sets of levels separated by 
integer multiples of 2tv/R with distinctive multiplicities. If these are absent, then the defect 
D cannot be conformal. 

For the two special cases of a topological defect and a factorised defect the spectrum 
takes more particular forms. 

If the defect D is a topological defect then the partition function is simply that of the 
original model on a strip of width R with two conformal boundary conditions, so that 

Z(B,B');D topological = y £ j rn i Xi,c{Vq) ■ (5.20) 

i 

In this case, the spectrum falls into sets of levels separated by integer multiples of tt/R, 
which is half the spacing of the energy levels in the general case. 



If the defect D is in fact a factorised defect of the form (2.12) then the partition function 
takes the special form 

Z (B,B');D facto „ sed = Yl ( n ^Nha N B'b) Xi ,MXj M • (5-21) 

In this case, the spectrum falls into sets of levels separated by integer multiples of 2tt/R 
as in the general case but with typically much higher degeneracies than those of a general 
conformal defect. 

A further check can be obtained by moving the defect across the strip to different angles 6. 
For topological defects the spectrum should be independent of 9 while for purely factorising 
defects F = ^ a b \ \a)) the partition function is 

Z (B,B>y,D factoris Jq;v)= E ( n ohNh a Ni, b )xi,c(^)Xj,M^)^ d = V*, (5.22) 

a,6,i,j 

which lets us identify the coefficients n a b exactly. 
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6 TCSA results 



We consider two different systems: firstly the Ising model, as a check on our method and 



confirming the results of | FFN |; secondly, we consider minimal models with p large for which 
we take p = 10 as a typical example. Throughout this section we use the dimensionless 
couplings appropriate to the TCSA method. 

6.1 The critical Ising model 



Using the defect TCSA method we can check the predictions of [FFN] as well as the accuracy 



of the TCSA method. We recall that fFFNf predict flows of the form X a -> iV(5vr/4 - a) 
where 

tan(a) = K r /ni , (6-1) 
As an example, we calculate the TCSA spectrum for a = 3ir/8 and compare it with the 



partition function at ip$ = 57r/4 — 3tt/8 = 7tt/8 given in (3.7) 

1\ _J^L_J_/ I 5 9 9 

Zd(0) + ,d(77t/8)W = 1 128 24 ( 1 + Q i + <? + <? 4 + 2<f + 2g4 

n 13 15 , 17 9 19 K x 

+ 3q 6 + 3q~ + q~ + 5g 4 + 5q~ + qz + q~ + 7<f + . . . ) . 
In figure 7(a) the normalised energy differences 2{Ei — Eq)/(E2 — Eq) are shown against the 



logarithm of the dimensionless coupling strength k of the perturbation. With this normali- 
sation, the indicator of an IR fixed point is the existence of a state with normalised gap 4. 
The only such point visible is that with log(ft) ~ 0.3. At this point, the energy levels visible 
have rearranged themselves into four distinct sets and both this arrangement and the values 
of the energy gaps agree with the partition function fl6.2|) . We conclude that the endpoint of 
the flow starting at the angle a = 3n/8 in the (m, K r ) plane is indeed the conformal defect 
,0(571-/4 - 37r/8) = D(7tt/8). 

The IR fixed point is located at a finite value of the bare coupling constant and at 
approximately the same distance from the origin in every direction. This means we can view 
all the IR fixed points together by plotting the energy differences ^{Ei — Eq) against the 
angle a at this fixed value of k. Expanding the partition function (|3~7|) gives 



Z D [0) + ,D{ Vo ){q\ \) = (1 + q 2 ~ 2x + q + q 3 ~ 2 * + 2q 2 + 2q 4 ~ 2x + q 2+2x + ...), (6.3) 

where x = (fo/ir. This means that the energy gaps are straight lines as functions of ipo or 
a and this is exactly what can be seen in the TCSA plot (figure |7(b)| ). The fact that these 



lines are straight also implies that k\/k,2 is a constant along the finite-size scaling flow. This 
is not expected to be the case in general, but the extra symmetry in the Ising model allows 
this to happen. 

This result provides a very strong indication that our identification of the RG flow fixed 



points in terms of Dirichlet-type defects is correct, confirming both the calculations of [ FFN ] 
and the usefulness of the TCSA method. 

6.2 Minimal models with p > 3 

For the minimal models M(p,p + 1) with p > 3 the flow picture is different from that for 
the critical Ising model, as is to be expected from the RG analysis. The proposed landscape 
of flows can be seen in figure [|. In the following we illustrate the exploration of this picture 
using the M(10, 11) model with TCSA. 
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(a) TCSA results in M(3, 4) for a = 3tt/8 plotted against (b) The spectrum of M(3, 4) at log(fi) = 0.3 as a 
log(ft) for N = 24 and 762 states. The degeneracies of the function of a for for it/4 < a < 57r/4 with N = 24 
energy levels at the IR fixed point are grouped according and 762 states, 
to the representation of the c = 1 folded model into which 
they fall. 



Figure 7: IR fixed points in the critical Ising model 



N 


10 


12 


14 


16 


18 


20 


22 


l°g( K *measured) 
lo g( K *thcory) 


-1.51 
-1.46 


-1.46 
-1.43 


-1.43 
-1.40 


-1.40 
-1.38 


-1.37 
-1.35 


-1.35 
-1.34 


-1.33 
-1.32 



Table 1: The approximate position of the perturbative fixed point in M(10, 11) in the purely 
chiral direction as a function of TCSA cut-off N 



6.2.1 Chiral perturbations 

The simplest cases are the purely chiral perturbations. The spectra for these models are 
both formally and numerically identical to the corresponding boundary perturbations as is 



explained in section 2.4. From section [37l| we see that the IR fixed points in the positive and 
negative directions are the (2, 1) and the (1, 1) defects, respectively. These defect flows are 
shown in figure || where the normalised energy differences are plotted against the logarithm of 
the coupling strength. In the positive direction the cut-off effects drive the flow past the first 
fixed point to a second one, realising the reversed version of the (1,3) — > (2, 1) flow. Since 
the spectra are numerically identical to those in the boundary case, this effect is identical to 



that observed for boundary flows in [Fc]. 



In the perturbative direction we can estimate the location of the TCSA fixed point as 
the value of k for which the fourth and fifth energy levels cross and this depends on N as 
given in table l[ As can be seen it is very close to the approximate value predicted by the 



analysis in [Fe| which gives k* ~ 0.1526A ra , once the different normalisation of the TCSA 



perturbing field is taken into account. 
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(a) positive direction: (1,2) — ► (2, 1) <■ — (1,3) (b) negative directional, 2) —+ (1, 1) 

Figure 8: M(10,ll): chiral flows. We show the normalised energy levels plotted against 
log(K) with N = 22, 1794 states. 

6.2.2 The third quadrant: k\ < 0, K r < 

Next we turn to the third quadrant of the plane of the couplings in figure |]. We have searched 
thoroughly for fixed points in the domain of convergence of TCSA and found a single fixed 
point F, which is in the diagonal direction ki = K r < 0. 



As discussed in section 5.1, the TCSA Hamiltonian for block diagonal and the 

flows for the two sectors are shown in figure ||. We find that for large k the spectrum 
approaches an IR fixed point for which the eigenstates in each sector separately fit into 
representations of the folded model. 

Combining the two sectors is a non-trivial task, because they have different effective 
cut-offs. This means that they have different values for the non-universal boundary energy 
and different re-scalings of the couplings k and r so that the two sectors cannot be directly 
compared. The usual perturbative renormalisation based on minimal subtraction or analytic 
continuation cannot be implemented within the TCSA framework, but it is possible to define 
a prescription for the relative shift of the sectors. Going from a TCSA cut N = 2n to 
N = 2n + 1 the even energy levels are not affected, because the new states in the Hilbert 
space belong to the odd sector. Similarly, raising the cut from N = 2n + 1 to N = 2n + 2 
the odd levels do not change. One can define, say, the values of the even eigenvalues at the 
odd cut N = 2n + 1 as 

ef en {2n + 1) = 1 (ef en (2n) + e f en (2n + 2)) . (6.4) 

The justification of this prescription is provided by its application to the fixed point F where 
it yields the spectrum shown in figure ^. It corresponds to a partition function of the form 
(|5.22D with rj = 1/2, namely 



1 9 

zK<z;«) = Z>,i(g) 2 . (6.5) 



r=l 
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(a) Even sector 



(b) Odd sector 



Figure 9: M(10, 11): diagonal flow in the negative direction. In the main figure, the scaled 
energy gaps for both sectors are shown, the odd sector with TCSA cut 19 and the even sector 
given by the combination of cuts 18 and 20 as detailed in the text. The small sub-figures 
show the even and odd sectors separately. 

although we could not identify all the components as some of them enter only at high level. 
This implies that 

the IR fixed point in the negative diagonal direction is the factorising defect 

p-i 

F = ^||r,l»((r,l||. 

r=l 

We found the same result for other minimal models and we checked also that the dependence 
of the spectrum on the position of the defect is as it is expected and given by ( |5.22j ). Note 
that the defect F is of the form of (|3.5| ) with s = s' = 1. 
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6.2.3 The second and fourth quadrants: kj and k t of opposite signs 



For different signs of the couplings k\ and K r we found strong indications of a fixed point 
far from the perturbative region which appears to be a deformation of the (//) fixed point 
in the Ising model; we call this conformal defect (//)'• However, since this is far from the 
perturbative region, there is no necessity that this be reachable by a finite-size scaling flow 
from the defect D, and indeed we believe that all RG flows starting at D in this quadrant 
end up at the identity defect /. This is the result predicted by the simple model in section [?] 
and to test it we explored the quadrant by making "round plots" similar to the Ising model, 
that is plots of the spectrum as functions of a at various fixed k for the Hamiltonian (5.6). In 



figure |10(a)| the weight differences are plotted against the angle a, defined in the same way as 
for the Ising model in ( |6.lD , with log(fv) interpolating between the approximate fixed point 
values log(/c*)= — 1.4 at 6 = tt/2 and log(K*) = —0.4 at 6 = it. The plot does not feature 
any fixed points but looks like a flow interpolating the chiral fixed points (2,1) — ► (1,1) 
located on the vertical and horizontal axes. This picture can be compared with an actual 
flow starting from the defect (2, 1). As can be seen in table fl4.9|) , in first order perturbation 



theory the field <fi transforms into the relevant field 



(13)(33) 



along the flow (1,2) (2,1) 



triggered by 4>. It is possible to study the flow starting from the defect (2, 1) generated by the 
non-chiral perturbation <^(X3)(33) in TCSA; some details of the method are given in Appendix 
A3. For negative values of the perturbation the flow is shown in figure 10(b) . It looks very 
much the same as the round plot, as shown in the schematic figure 111 




2.0 2.5 

(a) round plot 




0.2 0.4 0.6 0.8 1.0 

(b) (2, 1) - 13 ,33 HOW 



Figure 10: M(10, 11): comparison between the (2, 1) — ^1333 flow and the round plot. The 
absence of states of scaled energy 2 and 4 means that there is no indication of a conformal 
fixed point between the points a = n/2 and a = ir. The different slopes of the lines is a 
result of the different renormalisations of the function r(«) in the two systems. 

This tells us that: 

All RG flows starting from the defect D with Ki and K r of opposite signs end at the 
identity defect I, and the boundary of this space of flows is an RG flow from D(2,i) to I. 



6.2.4 The first quadrant k/ > 0, K r > 

Based on perturbation theory we expect a non-topological fixed point C in the diagonal 
direction K[ = K r > 0. We can apply the tools used in the previous cases for exploring the 
flows in the first quadrant. In figure 11 the normalised energy differences are plotted against 
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the logarithm of the coupling strength in the positive diagonal direction for the even sector. 
We find multiple degeneracies at log(re) ~ —1.2 indicating the presence of a conformal fixed 
point. One feature of this fixed point characteristic of a perturbative fixed point is that the 
energy levels stay close to their UV values, so that lines appear to split apart initially then 
to rejoin and intersect with multiple degeneracies at re*. 





-3 -2-10 1 

— r i/ 




(a) 4th group of lines (b) 5th group of lines 



Figure 11: M(10, 11): diagonal flow in the positive direction, even sector, scaled energy gaps 
vs log (re) for a strip of width tt truncated at level 20; the boxed regions are shown enlarged in 
the two smaller figures. The energy levels at the fixed point can grouped into representations 
of the folded model and are labelled m l where m is the degeneracy and I is the label of the 
representation. 

As shown in the magnified plots of the regions around the fixed point, the states approx- 
imately group into representations of the folded model with the correct multiplicities. An 
analogous pattern of line crossings is also seen in the odd sector. 

From the perturbative picture we also expect all flows starting with re/ > 0, re r > to 
flow to C and the boundary of this space of flows to be be a pair of flows starting from the 
two -D(2,i) defects generated by 4>[i^z) f° r a > ^/^ and by (/>(S3;13) f° r a < n /^- We have 
again checked this by comparing "roundplots" in the two-parameter space at fixed re and 
varying a with the one-parameter TCSA flows -D(2,i) + 0(13,33) and -D(2,i) + 0(33,13) > finding 
again qualitative agreement and support for the picture shown in figure [T^b. 
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D 



(2,1) 




(a) The conjectured exact flow from D(2,i) to I = 
and the set of flows starting from ^(1,2) f° r 
various values of a in the second quadrant. 



D (2,l) + <£i3,33 




(2,1) 



^33,13 



(b) The conjectured exact flow from D(2,i) to I = 
and the set of flows starting from -D(i,2) for 
various values of a in the first quadrant. 



Figure 12: The conjectured flows in the first and second quadrants 



As the level p of the model M(p,p+1) is reduced, the TCSA results become less convincing 
as the new conformal fixed point moves further away from the UV fixed point. The last model 
for which there is convincing evidence of the new defect C is M(5, 6); in model M(4, 5) the 
triple line intersection seen in figure pd| a is no longer seen in TCSA at the trunctation levels 
we have been able to achieve. It is possible that increasing the truncation level would show 
signs of a conformal fixed point, or the position of the conformal fixed point may have moved 
a large distance; at the moment we cannot tell. 



7 Conclusion 

For p > 3, the perturbative analysis carried out in section || and the numerical work in 



section g suggests the pattern of flows shown in figure |l3|b. We conclude that for p 3> 3 
there is at least one IR fixed point which is a conformal defect that is neither topological nor 
just a sum over conformal boundary conditions, namely the one in the K r = K[ > direction. 

The Ising model M(3, 4) behaves differently form the other minimal models and is treated 
separately in section |6.1| , with the pattern of flows |l~3|a. For p = 3, the fixed points that have 
been identified for large p all still remain, but are now just particular points in a continuum. 

For p small but greater than 3, the situation is not so clear-cut. We are still not sure 
exactly how the fixed-points approach the c = 1/2 continuum, which in part is due to the 
numerical inaccuracies of the TCSA method. 

The one-loop beta functions as discussed in section ||| predict the perturbative fixed points 
that we find with the TCSA analysis. It is tempting to try to extend this analysis to higher 
loop order to see if they confirm the RG flow structure we suggest in this paper. As a 
starting point, we consider the Ising model and then the changes that are present in the 
higher minimal models. Given the symmetries of the Ising model, the absence of three-point 
couplings and the continuum of fixed points, the RG flows away from the D defect can be 
modelled by the beta functions 
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(//)' 




v'(+-)u(-+) = C 



e = D' 



(++)u(— ) 



(a) The finite-size scaling flows in the Ising model 
showing the ring of conformal defects as a function 
of a with the identification of the factorised and topo- 
logical defects. 




. (//)' 



(b) The conjectured pattern of flows out of the 
defect fl( li2 j showing the 6 IR fixed points for 
M(p,p + 1) with p large 



Figure 13: Finite-size scaling flows 



where d is a scheme-dependent constant. This has the required 0(2) symmetry and a ring 
of fixed points at nf + = l/(2d). The minimal models M PtP+ i with p > 3 have non-zero 
three point couplings and so the simplest change to these beta functions to incorporate these 
is 

(:)-(:)-<:D-^-K:)- ™ 

This breaks the degeneracy and (for dy > 0) has six non-trivial fixed points with the same 
pattern of flows that we have found, which is additional evidence in favour of the pattern we 
propose. 

There are also various physical quantities which we would like to calculate using the 
TCSA method which for various reasons have proved intractable, such as the expectation 
values of T and T (which would give a better test of a conformal defect than the examination 
of the spectrum) or the reflection and tranmission coefficients defined in flQRW ]. 

Finally, it appears likely that a qualitative description of the space of defects and defect 
flows can be found from the microscopic RSOS model description of the minimal models, 
something we plan to address in the future. 
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A Appendix 

A.l Defect operators commuting with -D( r ,i) 

Let 8(s,s') be 1/2 if s = s' = and 1 otherwise, as in the definition of F s \ s ' in (3^). To 
see that F s \ s commutes with the topological defects one simply computes 

p-i 



D (r,l)F 



s ,s 



E N {r})(tl) D ^)U,s)){{l,s'\\D {ttl) -5(8,8') 



t,u=l 



(A-l) 



^ 1, *»« 1, s' \\D (uA) D {r>1) ■ 8(s, a') = F s ^ s ' D 



M) 



u=l 



To establish the decomposition (3.5) of a conformal defect that commutes with -D(2,i) 

(and hence with all Di r W) it is enough to show that we can find constants in ( |3,6[ ) such 
that 

(C,||o))((6||) = (F,||o))((6||) for all a, 6 € 1 P . (A.2) 

Because of the identity ||r, s)) = ||p — r, p + 1 — s)) it suffices to verify this for factorising 
defects of the form || x, v)){{y, v' || with v < 2i- and t/ < 
For the left hand side of (IA.2T) we find 



(C7, «' ||) = (CT> (y ,i), || a?, w» « 1, t/ ||) = (C,D (y)1) \\x,v))((l,v' \\) 

p-i (A.3) 

a=l 

To reproduce this result with (F, \\ x, v )) ((y, v' ||) we will distinguish four cases, namely 
whether v and i>' are less or equal to , respectively. Of course, the three cases where 
either v or v' are equal to ^rp- can occur only for p odd. We will abbreviate 

(A.4) 

For ma s we make the ansatz 

m^/iMO-^lksMMl) , (A.5) 

where s') = ^ if at least one of s, s' is equal to q, and fi(s, s') = 1 otherwise. The rria s 
have the property that for all s, s', 

mf = mf: a and = m s Jl a . (A.6) 

The first equality is an immediate consequence of || a, q)) = \\p — a, q)), while for the second 
one we need to commute the defect operators -D( p ,i) and C and use ((p, q\\ = ((l,g||, 

m^ a = l(C,D {Pil) \\a,8))((l,q\\) = \{CD M , \\a, s))((l, q\\) 

(A.7) 

= i(C,||a, S ))«l, g ||£> (P)1) ) = §(C,||a,*»«l,g||) = . 
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From their definition it is easy to check that the factorising defects F s \ s ' obey the identities 

FP +i-s\p+l-s> = F s\s> ; D (ajl) F p+1 - s \ s ' = D ip _ aA) F s ^ . (A.8) 



This allows us to restrict the range of the s and s' summation in (|3.6| ) to s, s' < q. We will 
now establish the equality (A. 2) in the four cases. 

i) v, v' < q: Then 

L<?J p-i 

(F,\\x,v))((y,v'\\) = E E < |s '^^(^)%i)l|l^))((l^'ll%i) ) lk^))((y^ , |l) 

s,s'=l a,t=l 

LqJ p-i 



s,s'=l a,t,u=l 
p-1 



(as,!) 



= ^(C,||a,,))((l,,'||)iV (ail)(j/il) , 

a 

(A.9) 

where in the last step we used that since s,v < q, \\u,s)) = \\x,v)) if and only if u = x and 
s = v, and similarly for ((t, s' || and ((y, v' ||. Since v, v' < q this implies also 5(s, s') = 1. The 
above expression is equal to QA.3| ) by the symmetries of the minimal model fusion rules. 

ii) v < q, v ' = q: Then 

LsJ p-i 

{F,\\x,v))((y,q\\)= £ £ m^ s 'iV ( Jjg/^ «') (||«, a » , || z, «» ((y,q\\) 

s,s'=l a,t,u=l 

a=l a=l a=l 

(A.10) 

where we used that N^ a i^_ y y = ^( p -a i)( y i ) • Since mp_ a = mj', the factor of \ in ( |A.5| ) 
ensures that both sums combine to give (A. 3). 

hi) v = q, v ' < q: This case works along the same lines. 

iv) v = q,v' = q: Then 

LsJ p-i 

[F,\\x,q))((y,q\\) = E E ^iV ( J^ ) tf( a>fl ')(||«, a >><<t lfl '||,||x, g >><<y 1 9||) 

s,s'=l a,t,u=l 

P-1 

= y mt ■ \ • (n, §\ + N, + N, Ifr*'? + iV, ,t' x '\) (A.11) 

^ a 2 V (a,l)(j/,l) (a,l)(p-y,l) (a,l)(j/,l) (o,l)(p-j/,l)/ 

a=l 

a=l 

where we used that N, if, ,., = -W, if/ The rest of the argument is as in ii). 
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A. 2 The renormalisation group and finite-size scaling relations in TCSA 

In this section we derive the form of the finite-size scaling flow in TCSA given in equations 
( |5.15 )-( 5.17l ). We start from the three assumptions in section |oT^ : if K{(N,t) are the TCSA 
coupling constants along a finite-size scaling flow parametrised by t at cut-off N, then we 
assume 

Ki(oo,t) = e yt K°i , (A. 12) 

~ N M (N,t) = PMN,t);N), (A.13) 
&(k;JV) = N v ji(K,N~ v ) . (A.14) 

Firstly we remove the expected t dependence from Ki ( Ne l , t) and define 

a\N, t) = e- yt Ki{Ne\ t) . (A.15) 
This satisfies the differential equation 

- n ^\n, t) = Ny 7i (N-y<j) , (A.i6) 

with initial conditions 

cr f (oo,£) = k° . (A.17) 

Since both the initial conditions and the differential equation are independent of t, the 
solution is also independent of t and we find that 

ai(N,t) = fi(K° i ;N) . (A.18) 



Substituting this into ( A.15 ) we find our first result, 

m{Ne\t) = e^fi^-N) = e^(A,0) . (A.19) 
Rescaling N in this equation gives the second result 

Ki(N, t) = e yt Ki(Ne- 1 , 0) . (A.20) 
Finally, differentiating this with respect to t at fixed level N$ gives 

^i(N ,t) = ye y t Kl (N e- t ,0)-ey t N e- t ^(N e- t ,0) 

= y Ki (N ,t) +e^(« i (iVbe- t ,0);JVoe-*) 

= y Ki (N , t) + NVjiiN-VeVtKiiNoe-t, 0)) 

= y Kt (N , t) + k^Ki{N Q e-\ 0); N Q ) 

= yK i (N ,t)+f3 i (K i (N ,ty,N ) 

= 0i(K(No,ty,No) . (A.21) 
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A. 3 Some details of the TCSA algorithm 

In order to calculate the matrix elements of the perturbations one needs to know the func- 
tional form of the 3pt function of primary fields on the upper half plane. If the perturbing 
operator is chiral it is simply 

(A\<f> h (z)\B) = Cz hA ~ h - hB . (A.22) 

If the perturbation is non-chiral then the correlation function is more complicated, but if 
the states A or B belong to the representation (1, 2) or (2, 1) the null vector equations can 
be used to deduce the form of the correlator. 
The equation in the model M(p, q) for B is 

(A\cj ) (z,z)(L 2 _ 1 -tL^ 2 )\B) = (A.23) 

where t = p/q if B = 4>i,2 and t = q/p if B = 02,1- Commuting the L n 's past cj>, changing to 
polar coordinates (r, 9) and using 

(h A -h B )(A\<j> hrh (r,e)\B) = {A\[L Q ,<j> hJl {r,9)\\B) = (h + h + r^) (A\</>(r, 9)\B) (A.24) 



one arrives at the second order differential equation 



[sin 2 6 d 2 + (1 + A - t) sin 29 d e + 

A(A + 1) cos 2 9 — A sin 2 9 - t{A cos 29 + h e~ 2ie + h e 2ie )] {A\(j) h - h (r, 9)\B) = (A.25) 

where A = h^ — h — h— hs- 

The null vector equation for A reads 

((L 2 _ x - tL_ 2 )A\cl>(z, z)\B) = . (A.26) 

With the same steps this can be converted to the equation 

sin 2 9 dg + 2sin9 ((1 + 2h + 2h - t - A)cos9 + 2i(h -h)sin9) d e + 

e 2ie ( 4/l 2 _ h ( 3t + 2A - 2)) + e~ 2W {Ah 2 - h(3t + 2A- 2)) + 

i A(A + 2t-2) cos(20) + \{Ah- A)(4h - A)] {A\4> hJl (r : 9)\B) = . (A.27) 



If it is not true that = kg and h = h at the same time the equations (A.25) and 
( [A. 27 ) can be combined to obtain a first order equation. In case h^ = Kb and h one 
gets 

[d B + (1 + h + h - 2t) cot 9 + i{h- h)] {A\<f> hti (r, 9)\B) = (A.28) 
having the solution 

{A\<f> hj - h (r,9)\B) = Cr~ h -~ h e- i{h -~ h)e sm9- 1+2t - h -~ h = (A.29) 
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The matrix elements of the perturbation in the basis ( [5.3; ) can be calculated using the 
relations 

(L^ n A\cf>(z,z)\B) = ]T (V 1 ) (z n ~ k (A\L k ^z,z)\B) + z n - k {A\L k <p(z,z)\B)) + 
k=-i ^ ' 

(A\<t>(z,z)\L n B), (A.31) 

OO , I 1 \ 

<A|«Kz,*)|£-„B> = - ^ ("JV" (z-«- fc (^|L fc 0(z,z)|i?) + z- n - k (A\L k <l>(z,z)\B)) + 
k=-i V + / 

(L n A|0(z,*)|fl>. (A.32) 

For cj) left chiral the terms containing z are not present. Applying these relations iteratively 
every matrix element becomes a sum of correlators of the form 

(AlL^LL^z^^B) A,B primary. (A.33) 

The action of L_i's translates to derivatives of ( |A.22 ) or ( A.3C| ) with respect to z and z. 

A. 4 Position invariance of the spectrum for a chirally perturbed defect 



Let H{9) be the Hamiltonian in (5.2) for the perturbation by a single chiral field (j>h,o, 

H{6) = H + Hj(9) where H = ^ (l - ^) , Hj(6) = X l (1)* e lhe M (e ie ) . (A.34) 
We will show that 

H(0) = e ieLo H{0)e~ i9Lo (A.35) 

by verifying it on matrix elements. Clearly, e idL ° Hoe~ teL ° = Hq, so that it is enough to show 
Hi(6) = e tSL ° Hj(0)e~ teLo . On the one hand, from ( |5.8| ) we have, in the notation of that 
section, 

(v l \H I (e)\v J ) = A, e^c^-h-A,) ^ (A36) 

where the constant C is given by C = (v i\4>h(l)\vj) . On the other hand, 

(vile^Hj^e-^lvj) = e i6 ^- A ^(v t |#/(0)K) = e^-^A, (~f ( Vi \^ (l)\v,) 

(A.37) 

These two expressions coincide. We see that H{9) and H(0) are related by a similarity 
transformation, and hence have the same spectrum. For the application to TCSA it is 
important to note that the operator e i8L ° commutes with the projection Pjy to the truncated 
Hilbert space, so that also 

P N H(6)P N = e ieLo P N H(0)P N e~ i6Lo . (A.38) 

Thus even in TCSA, the truncated Hamiltonians for different positions of the chirally 
perturbed defect on the strip have exactly the same spectrum. For a perturbation by an 
anti-chiral field 0o,fr> the argument is the same, except that the similarity transformation in 
this case is H(0) = e~ i6Lo H(0)e ieL ° . 
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